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<^: 

, Abstract 

M , 

C^ , We estimate the couplings in the Heavy Hadron Chiral Theory (HHCT) la- 

^^ ■ grangian from the QCD sum rules in an external axial field. Stability of the sum 

C"*"- , rules at moderate values of the Borel parameter is poor that probably signals 

^"^ ■ a slow convergence of the OPE series. At large values of the Borel parameter 

, , they stabilize, and yield the couplings much lower than the constituent quark 

^ ' model expectations. This region is not trustworthy for baryons, but in the me- 

f^ — , son case only unexpectedly large contributions of a few lowest excitations could 

^O ' invalidate the prediction gi « 0.2. 

(N ■ 

o: 

^ ' 1 Introduction 

0\ ; 

i-G I It is well known that the QCD lagrangian with ni massless flavours has the SU{ni)L x 

I ■ SU{ni)ji symmetry spontaneously broken to SU{ni)v giving the {nf — l)-plet of 

Mh. pseudoscalar massless Goldstone mesons (pions) vr* (ttJ = 0). Their interactions at 

(— I \ low momenta are described by the chiral lagrangian (see e. g. |jl|) 

'O: L^ = ^Tra^S+a^S + ..., S = exp— , (1) 

H ; 

5t , where the pion constant / « 132MeV is defined by 

(e*- is the pion flavour wave function), and dots mean terms with more derivatives. 
Light quark masses can be included perturbatively, and lead to extra terms in (H). 
SU{ni)L X SU{ni)fi transformations act as S — > LSi?"*". Let's define ^ = expivr//, 
S = ^^; it transforms as ^ — > LSJJ~^ = US,R'^ where U is a SU{ni) matrix depending 
on 7r(x). The vector v^ = \{£,'^d^£, + id^^~^) and the axial vector o^ = K^+S^^ — 
^9^^"*") transform as v^ — > C/(f^ + d^)U^ , a^ — > Ua^U^. There is a freedom in 
transformation laws of matter fields such as V* because it is always possible to 
multiply them by a matrix depending on vr. The only requirement is the correct 
transformation with respect to SU{ni)y {L = R). It is convenient to choose tp -^ 
Uip. Then the covariant derivative D^ = d^ + v^ transforms as D^ip — > UD^ip. 
Covariant derivatives of tensors with more flavour indices are deflned similarly. 
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Hadrons with a heavy quark are now successively investigated in the framework 
of the Heavy Quark Effective Theory (HQET) Q (for review and references see [^ ) . 
To the leading order in l/rri, the heavy quark spin does not interact and can be 
rotated or switched off at all (spin- flavour and superflavour symmetry). The Qq 
mesons with a spinless heavy quark form the j = ^ n;-plet ip'^ . The Qqq baryons 
can have j^ = 0"*" or l"*" giving the scalar flavour-antisymmetric ni{ni — l)/2-plet 
A* and the vector flavour-symmetric ni{ni + l)/2-plet S . Switching the heavy 
quark spin on gives degenerate 0^ and 1^ B and B* mesons, ^ ^ baryons, and 
degenerate ^ and | S and S* baryons. 

Interaction of these ground-state heavy hadrons with soft pions is described by 
the Heavy Hadron Chiral Theory @, 0]- Excited mesons were incorporated in [E], 
and electromagnetic interactions — in [0]; chiral loop effects were considered in ]8|, 5]. 
We start from the HHCT lagrangian with the heavy quark spin switched off: 

L = L^+ IPiiDoi^' + A*jiDoA'^ + ^*j- {iDo - A)S '^ (2) 

+5i?,a;.75V'^' + 2i52S*, • a) x S^''^ + 2^3 [A*,a] • S^''^ + S*, • a)A^'' 

where A is the S-A mass difference. The possibility of consideration of the SAvr 
interaction in HHCT relies on the fact that this difference is small compared to the 
chiral symmetry breaking scale though formally both of them are of the order of the 
characteristic hadron mass scale. The matrix elements of the axial current between 
heavy hadrons are easily obtained using PC AC: 

<M'\j]\M> = giu'ji^5u\ (3) 

<S'|/;|S> = 2^526^ X e*^ <A|4|S> = 2gse*i,e'\ 

where u^, e*-', e*-' are the M, A, S wave functions, and the nonrelativistic normal- 
ization of the states and wave functions is assumed. If we switch the heavy quark 
spin on, we obtain the usual HHCT lagrangian Q, |5|. 

The HHCT couplings gi should be in principle calculable in the underlying 
theory — HQET, but this is a difficult nonperturbative problem. Some experimental 
information is available only on gi. If we neglect l/rric corrections, then T{D*~^ —f 

n 71 
D TT^) = ^, and similarly for D^ir^ (with the extra 1/2). The experimental 



upper limit [^ on T{D*~^) combined with the branching ratios ||10[] give gi < 0.68. 
A combined analysis of D* pionic and radiative decays was performed in ; it gives 
g^ ~ 0.4-0.7. 

In the constituent quark model, gi is the axial charge g of the constituent light 
quark in the heavy meson. Moreover, following the folklore definition "constituent 
quark is B meson minus b quark", this is the most clear way to define g of the 
constituent quark. The baryonic couplings 52,3 sue also equal to g in this model. 
The most naive estimate is g ~ 1; the nucleon axial charge is gA = §5 in the 
constituent model, and in order to obtain gA = 1.25 we should assume g = 0.75. 

Sum rules [|ll| were successfully used to solve many nonperturbative problems 
in QCD and HQET. The currents with the quantum numbers of the ground state 
mesons and baryons with the heavy quark spin switched off are 

fM = Q*^-^Q\ fli = {q''^'CjW^)Q, A = {q'^'C^q'^)Q, (4) 

where Q is the spinless static quark field, C is the charge conjugation matrix, q^ 
means q transposed, {ij) and [ij] mean symmetrization and antisymmetrization. 



There are also currents jji^2 ^^'^ J*e2 ^i^h the additional 70. Correlator of the mesonic 



currents was investigated in |12, 13|, and of the baryonic ones — in [14|. The sum rules 
results are in a qualitative agreement with the constituent quark model: the massless 
quark propagator plus the quark condensate contribution simulate the constituent 
quark propagation well enough. 

The sum rules method was generalized to the case of a constant external field for 



calculation of such static characteristics of hadrons as the magnetic moments |15]. 
Sum rules in an external axial field were used |1^ for calculation of qa of light 
baryons. In the present work we use HQET sum rules in an axial field to calculate 



51,2,3- Sum rules for the D*Dtt coupling have been studied earlier ||171]; their result 



(in the terms used here) is gi = 0.2. We disagree with some results of this paper. 

2 Mesons 

We introduce the external axial field A^- {A] = 0) by adding the term 

Ai = j},Ai, (5) 

to the lagrangian. We are going to calculate correlators of the currents (Q) up to 
the terms linear in A (these terms are denoted by the subscript A) . The light quark 
propagator in the gauge x^A^{x) = gets the contribution 

The G^ term in Sa vanishes after the vacuum averaging; we are not going to calculate 
gluonic contributions beyond G^ and hence may omit this term.^ The axial field 
induces the quark condensate^ 

<r''ix)q^bp{0)>A = §^{[fA]+'-<qq>{[A),x]+6A^- 



2f2 



mjf 
36 



5x'^Ai-2Aj-xx) 75}°, (7) 






where mf ^ 0.2GeV^ is defined by [0]0 <0\dgG fj,,yju75u\'K^ > = imffp^, N 

In 



is the number of colours, Cp = -t^t"- We assumed p ■ A = where p ^ is the 



A-p 
The correlator of the meson currents (Q) has the A-term 



momentum of the field A\ in general these formulae should contain A^ = A 'fp. 



i<Tf{x)j^{0)>A = ^-^A]j,^^6{x)Ua{xo) (8) 

that depends only on A* and not A*q. The correlator possesses the usual dispersion 
representation at any A^ . The meson contribution at Aj = is p{u)) = F'^6{uj — e) 

*The first term in (n) was presented in [M, but we don't understand how the authors calculated 
G^ corrections without the second term and the above statement about the third one. 

^ These formulae were presented in |lq] but with a different second term in the first equation; 
their term does not satisfy the relation <{Dq)q> = that follows from the equation of motion. 

®the sum rule considered in the second paper of llq] yields the relation mf ~ mo/4. 



where e is the meson energy and <0|j*|M> = iFu'^ (the usual meson constant is 
/m = -n=)- Switching on the external field produces the energy shift e — > £—51757- A 

(sn = 2757'"' ~ ^2 ^^ ^^^ meson spin projection). Therefore the meson contribution 
is /3a(^) = giF'^5'{ijj—e)+c5{ijj—e) where the second term originates from the change 
of F^. Besides that there is a smooth continuum contribution p™'^*(a;). Thus we 
obtain 

liAiu:) = /^ + ^ + nr (^)- (9) 

[u) — e)'^ uj — e 



In other words, the lowest meson's contribution in both channels (Fig. |T|a) gives the 
double pole at u; = e; mixed lowest-higher state contributions (Fig. |l|b) give a single 
pole at a; = e plus a term with a spectral density in the continuum region after the 
partial fraction decomposition; higher states' contributions (Fig. ||c) have a spectral 
density in the continuum region only. 

We can also calculate the correlator using OPE. Gluons don't interact with the 
heavy quark in the fixed point gauge x^A^(x) = 0. The diagrams with non-cut light 
quark line (Fig. ga, b) vanish. The diagram with cut line (Fig. ^) gives 

rr f+\ f^ f-t ^<gg>^ 5 2,2\ /in\ 

^Ait) = ^ 1 :rH ^"^1* • (10) 



4 V 3/2 36 

Thus the appearance of gA_ of the constituent quark is entirely caused by interaction 
with the quark condensate. The correlator (|lO|) corresponds to the spectral density 
in the form of S{uj) and its derivatives. Equating these expressions at an imaginary 
time t = —i/E we obtain the sum rule 

'"^%c)e-./- = 4(1-^ + 4^). (U) 



E J 4 V 3PE 36^2 

At sufficiently high energies the hadronic spectral density is equal to the theoretical 
one, and it vanishes. Moreover, higher states' contributions are exponentially sup- 
pressed, therefore we neglect them. This sum rule should agree with the rric — *• 00 
limit of the sum rule for the D* -^ Dn coupling ||l^ . We disagree with the last term 



in [17], though it is not very important numerically. 



At large E the sum rule ( pT| ) leads to 

giF^ = \fh - ^<qq>. (12) 

Of course, higher states' contributions are not exponentially suppressed at large E. 
But the contribution of high energy excitations is given by the perturbative spectral 
density which is zero in our case. Hence only one or few lowest states contribute. 
Excited mesons are supposed to have F^ substantially smaller than the ground state 
meson, and may be are not very important. So we may hope that ( |l2|) gives us a 
reasonable estimate. We use -F^ = (240MeV)^ and e = 430MeV obtained from 
the sum rule without radiative corrections 1131 because we don't know radiative 



corrections to (IC) and feel that it would be inconsistent to take them into account 
in a part of sum rules (this would lead to large errors e. g. in the sum rule for the 
Isgur-Wise function). Of course, the existence of large radiative corrections ||l^ leads 
to an additional uncertainty in the result. If higher states' contributions to ( p^ are 
not important, we obtain gi ~ 0.2. 

We can multiply ( [Tl| ) by exp(e/£') and differentiate in E in order to exclude 
c. The result is shown in Fig. H. Of course, at large enough E > 2GeV dl^ 



is reproduced, and we get gi ~ 0.2. This is the only region in which the sum 
rule is stable. A similar sum rule with the finite rric was analyzed in |]l^] in the 
region E = 1.5-3GeV with the same result. But this requires an extremely slow 



convergence of the OPE in order to forbid using lower E. The OPE series (11) looks 
like 1 + 200MeV/ii^ + (170MeV/£^)^, so it seems that the expansion is apphcable 
at £" > SOOMeV. This is close to the lower bound of the applicability region of the 



ordinary sum rule |12|. The upper bound of its applicability region is about SOOMeV 
because the continuum contribution is too large at higher E. If we suppose that the 
applicability region of the sum rule (|l^) is approximately the same, then we obtain 
gi ~ 0.4-0.7. Stability of the sum rule in this region is poor. In this case we have 
to assume that the higher states' contribution is significant at ii^ > SOOMeV. 

3 Baryons 

Correlators of the SS and AS currents have the ^-terms 

i<Tjl{x)r^,y{Q)>A = Al,5p{x)IiA{x(,). (13) 

If we define <0|j^|A> = F\e^, <Q\j^^\Tj> = -Fse^, then the physical states' con- 
tributions to the correlators (Fig. |l|) are 

(w-es) oj - et, 

(w -eA)(w -es) "^-^A w-ee 

It is impossible to separate the g^ term from the mixed A-excited and S-excited 
contributions unambiguously. We can do it approximately if A = e^ — sa <C ^c— eA,s 
because in such a case partial fraction decomposition of the first term would give 
large contributions to ca,s ~ 1/A with the opposite signs while the natural scale 
of CA,s in (|T^) is ca,s ~ l/(£c — ^A,s)- This is not a defect of the sum rule but 
the uncertainty inherent to g^ which can be defined only when A is small compared 
to the chiral symmetry breaking scale. We choose to require c\ = cs; the choices 
CA = or cs = would be equally good. 

The baryonic correlators in the OPE framework are described by the diagrams in 
Fig. ^. The diagrams Fig. ^-c with the non-cut quark line interacting with the axial 
field vanish due to (y). We use the factorization approximation for the four-quark 
condensate in Fig. ^. In this approximation two diagonal correlators <jiji> and 
<J2J2> coincide in both ES and AE cases, as well as <jaiJS2> and <JA2Jt;i>- This 



is similar to the usual heavy baryon sum rules |14], and confirms the observation 
that F\i = F\2 and i<x;i = i^S2 within the factorization approximation to the sum 
rules. Only even-dimensional condensates contribute to the diagonal SS and the 
nondiagonal AS correlators: 






NttH^ 



l-^±^ ^^ + 



5 Cb\ mft^ TT^<qq>^f 



3 CfJ 12 QNf^ 



(15) 



where Cb/Cf = 1/(A^ — 1) (this term comes from the diagram Fig. ^). Only odd- 
dimensional condensates contribute to the nondiagonal SS and the diagonal AS 



correlators: 



2JV! <qq> 



,2 j?2+2 



2A^ 



(1 + l^mlt" - ^mlt") 



These correlators correspond to the spectral densities 



Pa{^) 






^2 + ( :i ^ -B \ mi 



Cf 



Pa{^) 



2N\ <qq>u 

3A^7r2 



(16) 



(17) 



(plus 8{ijj) and its derivatives). 

We use the standard continuum model /9™°*(a;) = p^^^°^'d{uj — £c)- Equating the 
OPE (O, [16|) and the spectral representation at i = —i/E, we obtain the sum rules 



2g2K' 
E 
4/2 



^+c]e-^-/^ 



E' 



vr^ 






4<gg> ^2 



f2g^F\Fh 



TT f^ / 771 ^771 \ 



A 



, A c 

th h - 

2S 2 



e--A/i? + g-e5:/i? 



(18) 



^^ 



7mf 7r^<gg>^ 

J2 (ec/^j + 7^^/o (ec/-Ej + ^g „2^4 



4<gg> ^2 



TT' 






Where /„(rr) = l-e-E:;=o^- 

Fig. ^a shows the results for (72,3 obtained from the diagonal SS and the nondiag- 

onal AS correlators (the first formulae in (|l^)). The values of eA,s and -Fa,s are taken 

from the diagonal sum rules |14] in the middle of the stability plato E ~ 300MeV. 

We use the same values of the continuum threshold as in [^: (5.6 it 0.5)k for S 

and (4.6 ± 0.5)k for A where k^ = —\<qq>, k = 260MeV. We put the continuum 

threshold in the sum rule for §2 (SS correlator) equal to that in the S sum rule; in 

the case of 53 (AS correlator) we put it equal to the average of the thresholds in the 

A and S channels. Of course, the effective continuum threshold for the ^d-term in the 

correlators does not necessary coincide with that in the absence of the external field; 

this assumption gives an additional uncertainty. The spectral density ( p^ ) behaves 

like ijp , and hence the continuum contribution to the sum rules ( |l8|) quickly grows: 

it is equal to the result at i? ~ 600MeV and is three times larger than the result at 

E K, 900MeV. If we assume that the uncertainty in the continuum contribution is, 

e. g., 30%, then we can't trust the sum rules at £■ > 900MeV. The lower bound of 

the applicability region is determined by the convergence of the OPE series for the 

correlators. It behaves hke 1 + {l9f)M(N /Ef + (245MeV/£^)^ for the SS correlator; 

in the AS case 140MeV enters instead of 190MeV. It seems that OPE should be 

applicable at E^ > 400MeV. Stability of the sum rules in this window is poor. We 

can only guess that §2 ~ 0.2-1, g^, ~ 0.1-1, and probably 53 < §2- 

Fig. P^ shows the results for 5^2, 3 obtained from the nondiagonal SS and the 

diagonal AS correlators (the second formulae in (1^)). The values of eA,s and 

-^A,E are taken from the nondiagonal sum rules [|l4| in the middle of the stability 

plato E ~ 650MeV. Again we use the continuum thresholds from [ij]: (5.6 ± 0.5) A; 



for E and (3.65 ± 0.5)A; for A. The spectral density (|17D behaves hke oj, and the 
continuum contribution grows not so quickly: it is equal to the result &t E k. IGeV 
and is three times larger than the result at £^ ~ l.GGeV. The OPE series behaves 
like 1 + (290MeV/^)2 [l - (ISOMeV/S)^] , and the apphcability region starts at a 
larger E. Stability of these sum rules is poor too. They tend to prefer somewhat 
larger values 52 ~ 53 ~ 0.5-1. 

In conclusion, poor stability of the sum rules at moderate values of the Borel 
parameter probably signals a slow convergence of the OPE series. At large values 
of the Borel parameter they stabilize, and yield the couplings 51,2,3 much lower 
than the constituent quark model expectations. This region is not trustworthy for 
baryons due to large continuum contributions. In the meson case, higher excited 
states are dual to the vanishing theoretical spectral density, and don't contribute. 
Only unexpectedly large contributions of a few lowest excitations could invalidate 
the prediction gi k. 0.2. This inconclusive picture is especially surprising because 



similar sum rules in an external axial field produce qa = 1-25 for the nucleon [16| in 
excellent agreement with experiment and the constituent quark model. 
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Figure 1: Physical states' contributions to the correlator 



Figure 2: Correlator of meson currents 



Figure 3: Sum rule for gi 



Figure 4: Correlator of baryon currents 



Figure 5: Sum rules for §2^3 (solid and dashed lines): a — diagonal EE and nondiag 
onal AS sum rules; b — nondiagonal SS and diagonal AS sum rules. 



